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We evaluate the exact QED 2+1 effective action for fermions in the presence of a family 
of static but spatially inhomogeneous magnetic field profiles. This exact result yields an 
all-orders derivative expansion of the effective action, and indicates that the derivative 
expansion is an asymptotic, rather than a convergent, expansion. 

The effective action is a fundamental tool for the study of quantum field theory. 
Using the proper-time technique, Schwingeitl showed that the QED effective action 
can be computed exactly for either a constant or a plane wave electromagnetic field. 
This result was later adapted to QED 2+1 with constant fields by Redlicha. For 
general electromagnetic fields the effective action cannot be computed exactly, so 
one must resort to some sort of perturbative expansion. A common approach is the 
derivative expansioncrQ in which the field is assumed to vary only very 'slowly'. A 
useful computational procedure for the derivative expansioroQ is provided by t 
supersymmetric path integral representation of the proper-time effective actio 
However, even first-order derivative expansion calculations of the effective action 
are rather cumbersome, and it is very unclear just how convergent (or otherwise) the 
expansion is. Thus, the physical significance of a first-order correction term must be 
examined with care. In this paper we seek a better nonperturbative understanding 
of the derivative expansion by considering a new exactly solvable model which has a 
spatially varying magnetic field. The exact integral representation for the effective 
action may then be expanded in an asymptotic series, and we find that the first few 
terms agree with (independent) derivative expansion computations. 

Specifically, we show that the QED 2+1 effective action for massive fermions can 
be computed exactly^ in the presence of time-independent but spatially varying 
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magnetic fields of the form: 



B(x,y) = £sech 2 (-) (1) 



Here, the magnetic field is uniform in the y-direction but is inhomogeneous in the 
x-direction. The length scale A characterizes the width of the profile in the re- 
direction. In the limit A — > oo the magnetic field in (^) tends to a uniform field of 
strength B. The constant B (the peak value of the magnetic field) sets a length 
scale \j\l~B known as the magnetic length, and the derivative expansion regime 
corresponds to A ^> 1/yB. In practice, the system will be considered in a box in 
the y-direction of size L. The total magnetic flux <E> is then finite: <E> = BXL/tt. We 
show that the effective action has a simple exact integral representation involving 
elementary functions for all values of mass to, width A, and strength B. 
The fermionic effective action is given by 

iS eS = Log Det (ilp - to) (2) 

where Ip = r ) u D v = 7" (d v + iA v ), and A v is a fixed classical gauge potential. The 
Minkowski gamma matrices 7^ satisfy the anticommutation relations {7 1/ ,7 <T } = 
—2g v<J = 2 diag(l, —1,-1), and in 2 + 1 dimensions the irreducible gamma matrices 
may be chosen as the 2x2 matrices: 7 = cr 3 , 7 1 = icr 1 , -f 2 = ia 2 , where the a % are 
the 2x2 Pauli matrices. Note that an alternative choice, 7 = — cr 3 , 7 1 = —icr 1 , 
7 2 = —ia 2 , can be viewed as changing the sign of the mass, m — » —to, in the Dirac 
operator. The system with effective action (|2j) is not parity invariant since a fermion 
mass term breaks parity in 2 + 1 dimensions a. However, in this paper we shall limit 
our attention to the parity-even part of the effective action (0): 

iS = ^Log Det [{ilp - to) (-ilp - m)] = ^Tr Log (jf + to 2 ) (3) 

The Tr Log may be defined by a "proper-time" integral representationS 

1 r°° ds 



iS 

z Jo 



(V s ^ 2 ) (4) 



From this expression we see that in order to compute the effective action in this 
proper-time form we need to compute the 'proper-time propagator' 

Tr(e- S ^ 2 ) (5) 



When the field strength is constant, this propagator may be computed exactly 
because the Fock-Schwinger gauge choice, A p = \x v F v>i , leads to (Ip) 2 being a 
quadratic operator. Thus the propagatorjn (JsJ) is essentially that of a harmonic 
oscillator, leading to the simple expressioncl 



S = 



A f°° ds 



8tt 3 / 2 /„ s 3 / 2 



LF|coth(|.F|s) - - 
s 



(6) 
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where \J-\ = ^JF IJ , V F^ V and A is an area factor associated with the degeneracy. 
Note that the 1/s term in the square parentheses in (|^) corresponds to the zero 
field case; thus, one actually computes the difference between the effective action 
with and without the background field. 

Ideally we would like to go beyond this constant field case and compute the 
effective action in an arbitrary background electromagnetic field. However, as this is 
not feasible, we consider a derivative expansionH^, which is a perturbative expansion 
about the constant field case. Consider, for example, a Taylor expansion of the 
gauge field 

= \x a F ail + X -x a x^d a F^ + ... (7) 

Here the derivative terms dF are in some sense smaller (see below) than the F terms. 
The operator Ip 1 now involves terms quartic in x, so the proper-time propagator (||) 
can be computed, at best, in perturbation theory. An efficient way to compute these 
perturbative contributions is to re-express the proper-time integraLrepresentation of 
the effective action as a SUSY quantum mechanical path integraMB. However, there 
is a proliferation of terms appearing in the perturbation expansion. For example, in 
the 2+1 dimensional case there are 6 different terms ^ in the first order correction 
to the constant field case. For definitencss we choose one such term by restricting 
to a static, but spatially inhomogeneous, magnetic field (no electric field). Then 
the first order correction to the effective action is □ 

For m -C B/m (i.e. for the fermion rest mass energy scale much smaller than the 
cyclotron energy scale set by the approximately uniform magnetic background) this 
correction simplifies to 

AS = -J— #C (5/2) / d^x^L. (9) 



V2(4tt) 2 16tt w 'J B 3 /2 

where Q{z) = XmLi n ~ z i s the Riemann zeta function and C(5/2) ~ 1.34. 

Computing higher order terms in the derivative expansion is messy and not 
particularly illuminating. Thus, to obtain a better understanding of the meaning 
of the expansion we seek an exactly solvable effective action with an inhomogeneous 
background field. To be specific, we consider static magnetic fields. Then the 
proper-time expression (|i|) becomes 

s 1 ^ ds 



4^ __ 



t Various inequivalent contractions of tjhe tensors F^ v , daFfi^, g^v and e a ^ 7 appear. In 3 + 1- 
dimensions there are many more termsQ due to the epsilon tensor e a py$. 
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after taking the trace over the energy. We must also sum over the two spin projec- 
tions ± arising in the Pauli operator m 2 — D 2 ± B, where B = Fn is the magnetic 
field strength. 

An alternative expression for the static effective action is obtained by performing 
the proper time integral in (|l^), yielding 

S=l^2 Tl V m2 ~D 2 ±B (11) 

which expresses the static effective action as a sum over the Dirac sea energies for 
single fermion statealJ (recall that the static effective action is minus the static 
effective energy). 

Yet another expression for the static effective action is the "resolvent form" 



r dE E2Ti 
4-7 2iri 



1 



i 2 - D 2 ±B-E 2 



(12) 



which reduces to the Dirac sea form ( |ll| ) after integration over an appropriate 
energy contour, and which may also be derived formally from the original expression 
iS = iTrLog(^5 2 + m 2 ) after an integration by parts in the energy trace. 

The three expressions (|l0|), (pi]), and (Jl^) are equivalent but one may be more 
useful than the others, depending on the form of the Pauli operator. For example, 
consider the uniform magnetic field case. The static vector potential can be chosen 
as A = (0, Bx), in which case the Pauli operator becomes 

m 2 - D 2 ± B = m 2 - (-^-\ + B 2 (^x - ^\ ±B (13) 

where k is a momentum variable corresponding to the free y direction. This Pauli 
operator has the form of a 1-dimcnsional harmonic oscillator with frequency oj = 2B, 
and thus has a purely discrete spectrum of levels known as Landau levels, each of 
which is infinitely degenerate. In fact, the degeneracy is given by the (infinite) total 
flux, which in a finite area A is BA/2tt. Thus the proper time propagator (||) is 

£ TV ^ e -(m 2 -D 2 ±B)s^ = BA R -m?s £ g e _ 2B(n+ l±I) s = BA ^ s 

± ± 71 = 

(14) 

which agrees with Redlich's result (JsJ) . The sum over Dirac energies in ( pr|) may 
also be performed explicitly because the spectrum is discrete: 



D 2 ± B = ^Jm 2 + 2B{n + 1/2) ±B 



47T 

± n=0 



B 3 / 2 A 
2V2tt 



, 1 m 2 \ 1 
2 C + 7777 ' 



2B V 2' 2BJ 



(15) 
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where ((z, a) = Y^—pin + a) 2 is the Hurwitz Zcta functionS. Using the integral 
representatio: 



r(z) 



C(z,a+1) + 



a 



/ dse~ 2as (2s) z " 1 coth s (16) 
Jo 



we see that the proper time expression (^) provides an integral representation of 
the Zeta function expression (jL5j). 

These two approaches [the proper-time method ( |Io|) and the Dirac sea method 
(|TT|)1 are the most direct, and most familiar, ways to compute the effective action 
with a constant B held. However, in the case of an inhomogeneous magnetic field 
considered below the resolvent method proves more useful. Thus we first review the 
resolvent approach to the constant field case. 

With the vector potential A = (0, Bx), the static effective action ( |l2| ) is 

Here is the resolvent Green's function for the one-dimensional Schrodinger 

operator in (p"3|). The diagonal resolvent may be constructed, as usual, as the prod- 
uct of two independent solutions u\(x) and U2{x) to the corresponding homogeneous 
differential equation, divided by their Wronskian: 

G(x,x) = ^p^l (18) 

The Wronskian W = u[(x)u2(x) — ui(x)u' 2 {x) is independent of x. For the harmonic 
oscillator operator in (Q) the independent solutions are parabolic cylinder functions 
D v (x), and so the trace of the diagonal resolvent gives (with a straightforward 
regularization) 

§ T ^i = ^(- ± ) 

where i/± = (E 2 - m 2 - B(l ± 1))/{2B), and ij}(z) = (d/dz)logT(z) is the Euler Psi 
functiorO. The Psi function has simple poles at the negative integers (correspond- 
ing to the familiar discrete spectrum of the HO system), and so the E integral in 
( |l7| ) may be performed to give the Zeta function expression in ([l5|) . 

To go beyond the constant field case, we choose a vector potential of the form 
A = (0, a(.x)), which corresponds to a magnetic field B{x) = a'(x) only depending 
on the x-coordinate. The Pauli operator for such a vector potential is 

m 2 -D 2 ± B(x) = m 2 - ( 4~ ± \a(x) - k] \ ( 4~ T [atx) - k]\ (20) 

\dx J \dx J 

which has the form of a one-dimensional"'' SUSY quantum mechanics operatorEl. 
This special form for the vector potential is chosen because for various particular 

■i-Note, however, that the super-potential, Qfe(x) = a(x) — k, is k dependent, where k is the mo- 
mentum for the free y-direction. This reflects the fact that the Pauli operator is really a two- 
dimensional operator. 
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choices of the superpotential, Qk(x) — a(x) — k, the SUSY Pauli operator has an 
exactly known spectrum. Given complete knowledge of the spectrum it should then 
be possible, at least in principle, to compute the effective action using one of the 
three expressions (|l0|), (11), or (p^). With Qk(x) = Bx — k we obtain the constant 



field case discussed above. With the superpotential choice 

Q fe (x) = BAtanh(^) - k (21) 
A 

the corresponding magnetic field profile is as in (|l]) , and the Pauli operator has an 
exactly solvable spectrurnciL 

However, there are two crucial differences from the constant field case. First, the 
spectrum depends explicitly on the momentum parameter k. Therefore, the trace 
over k does not simply lead to an overall flux factor [as was the case in (|lj) , (|l^) and 
fl9|)1, but must be truly integrated over. Second, the spectrum is no longer purely 
discrete - it has continuum states (both one-way and two-way scattering states) as 
well as discrete bound states. Moreover, the density of states depends explicitly on 
k, so that, for example, the actual number of bound states depends on k. Thus, an 
explicit sum over the spectrum as in ( [ll"|) is problematical, and it is nontrivial to 
find a simple representation of the proper time propagator in (|Io|). Nevertheless, it 
is possible to find the trace of the diagonal resolvent in ( |T^ ) . 

The corresponding homogeneous differential equation can be converted to hyper- 
geometric form with the substitution £ = [1 + tanh(x/\)]/2. The two independent 
solutions may be chosen as (here ± refers to the two spin projections in (|2p|)) 

M x ) = T(l - €)fF(ot + P T B\ 2 ,a + p + 1 ± BX 2 ; 1 + 2a; £) (22) 

u 2 {x) = e*(l - O F{a + (3 T BX 2 ,a + + 1 ± BX 2 ; 1 + 2/3; 1 - £) (23) 

Here, a and (3 are defined as the following roots (with positive real part): 

a=y a 2 _ k -E 2 {3=y a l-E 2 a k = V(k-BX) 2 +m 2 (24) 

With these conventions, u\ is regular at £ = (i.e. as x — > — oo), whereas u 2 is 
regular at £ = 1 (i.e. as x — * +oo). The Wronskian of these solutions is a constant: 

I / \ / \ , s s 2 r(l + 2a)r(l + 23) , . 

The trace of the diagonal resolvent is then 

J dxG { ^ k (x,x) =-L + [^(a + fiTBX 2 )+^(a + B+l±BX 2 )] (26) 

The spectrum contains a finite number of bound states, which arise from the (sim- 
ple) poles of the ^-functions in (|2^) . The same discrete spectrum may be obtained 
by solving the homogeneous equation for real E 2 and requiring normalizability of 
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the eigenfunctions0. The spectrum also contains a cut starting at a 2 and another 
cut starting at a 2 _ kl corresponding to the two barrier thresholds of the 'potential' 



m 2 + Ql± Q' in (M 



Remarkably, the remaining integral over k may be performedu, leaving an inte- 
gral form for the effective action: 

S =AEf^ /2 Z(z2 ,~f, „ T BX 2 ) + ^(z + 1 ± BX 2 )] (27) 
AttV J ZQ J z z - z l + (Am) 2 



where zq = Ay {BX) 2 + to 2 . This expression can be rewritten in terms of elementary 
integrals using the following representation of the ^-functioncil: 

1 f°° t dt 

^ X) = lnX -T X - 2 J (t 2 +X 2 )(e 2 ^-l) (28) 

After summing over both spin projections in (|27j), we find that the first two terms 
on the right hand side in (|2^) contribute B-independent terms to the effective action 
and therefore may be dropped. For the third term, the z-integral in ( p7j ) can be 
carried out exactly and we obtain the following finite integral representation for the 
effective action 

s _ L r°° dt / a (A 2 to 2 + v 2 ) ^ Ato - iv | e c \ ^ 

AttX 2 Jq e 27r * — 1 \ v Am + iv J 

where c.c. denotes the complex conjugate, and v 2 = t 2 + 2itBX 2 . Expression ( |29| ) 
gives an exact integral representation for the effective action with the inhomoge- 
neous magnetic background ([!]), just as (g) gives an exact integral representation 
for the effective action with a uniform magnetic background. 

In the limit of vanishing mass, to = 0, the only relevant dimensionless parameter 
is -BA 2 , which is the square of the ratio of the width parameter A to the magnetic 
length 1/yB. Since the uniform field case corresponds to A » l/VB, we expect 
the derivative expansion regime to correspond to the regime in which _BA 2 — ► oo. It 
is straightforward to develop an asymptotic expansion of the integral (^9|) for large 
BX 2 when m = 0, yielding 

= LX{Bf/ 2 ^ 1 T(j - 3/2)r(j" + 5/2) 

8^ {AkBX 2 )i r(j + i)r(j/2 + i/4)r(3/4 - + 1 ' 



LX(B) 3 / 



15 1 



C(3/2) - — C(5/2) 



(30) 



16tt ' BX 2 

where Q(z) is the Riemann Zeta functionQ. 

We now compare this expansion with the first two terms j^f the derivative ex- 
pansion (with to = 0) obtained using the proper-time methodtl: 



S = d A x 



C(-l/2)r D , M3/2 , 1 15 W C /oJV5(z)] 



[BOr)] 3 / 2 + -=-— -^C(5/2) 



2 



V2tt L WJ V2(4tt) 2 16tt sv 7 ; [B(:c)] 3 / 2 



(31) 
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The first term comes from (|15|) and the second from (g). For the magnetic field 
profile B(x,y) — Bsech 2 (x/X) in (|l]), the x and y integrals can be done (as ev- 
erything is static, the overall time scale has been absorbed into the definition 
of the energy trace throughout). Then, using the facts that J dxsech 3 x = tt/2, 
J dxsechxtanh 2 x — 7r/2, and £(3/2) = — 47r£(— 1/2), we see that these first two 
terms of the derivative expansion agree exactly with the first two terms of the ex- 
pansion in (|30|). Thus, it is indeed natural to identify the expansion in (30) as 
an all-orders derivative expansion. This expansion, however, is an asymptotic ex- 
pansion rather than a convergent one. Indeed, as is typical of asymptotic series, 
the magnitude of the expansion coefficients decreases at first, then bottoms out 
and eventually tends to infinity. Also note that the sign of the expansion coef- 
ficients follows the pattern +, — , — , — , +, +, — , — , +, +, ... . The eventual 
doublc-plus/double- minus pattern is reminiscent of Stirling's asymptotic series for 
the Gamma function. 

For nonzero mass m, there is another dimensionless parameter, B/m 2 , which is 
the ratio of the cyclotron energy to the rest mass energy. A double expansion of 
29b yields 



_ Lm 3 \ ^2, 

3=0 



(2B\ 2 y 



E 

k=l 



(2fc + j-l)! 



& 



2k+2j 



with B n the 



Ah 



Bernoulli number 



(2fc)! 
§0 



(2k + j-i)(2k+j-l) 



2B 

m 2 



2k+ 3 



(32) 



Each power in l/(B\ ) corresponds to a 



fixed order in the derivative expansion of the effective action. The zeroth order term 
agrees with the B/m 2 expansion of the exact constant B field answeaa, while the 
1/(5 A 2 ) term agrees with the B/m 2 expansion of the leading derivative expansion 
contributionO. 

I conclude by commenting briefly on some natural extensions of these results. 
First, these results could be extended to 3 + 1 dimensions; although there one 
must be careful to include the proper charge renormalization, which is present 
even in the constant field casefil. Second, it would be interesting to compute the 
finite temperature effective action for this QED2+1 system. If one computes the 
zero temperature effective action using the proper time expression ( p^| ) then the 
extension to nonzero temperature is immediate^: 



Sett 



40F 



± 



ds 



Airs 



Tr 



cxp - (m — D ± B 



(33) 



Here (3 is the inverse temperature, fi is the chemical potential and the fourth theta 
function is0 



h(u\T) = l + 2j2(-l) 



cos (2nu) 



(34) 



n=l 



In the zero chemical potential and zero temperature limit, one simply omits the 
theta factor in the proper-time integral, giving ([l^). However, since this theta 
factor is independent of the magnetic field, if one had computed the proper-time 
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propagator for this system (rather than using the resolvent technique described 
in this paper) then the exact finite temperature effective action would follow im- 
mediately. It would be interesting to investigate the properties of the derivative 
expansion at finite temperature using this exactly solvable model. 
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